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Abstract. The paper presents some results of qualitative analysis of Kirchhoff's differential 
equations describing motion of a rigid body in ideal fluid in Sokolov's case. The research 
methods are based on Lyapunov's classical results. Methods of computer algebra imple- 
mented in the computer algebra system (CAS) "Mathematica" were also used. Combination 
of these methods allowed us to obtain rather detailed information on qualitative properties 
for some classes of solutions of the equations. 
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1 Introduction 

' 

Let us consider the problem which may be of interest both for the development of the method 
. of investigation employed and for numerous applications. In the problem proposed for conside- 

ration the equations of Euler-Poisson's type are the model of the object under investigation. 

Euler-Poisson's differential equations describing the motion of a rigid body with one fixed 
point and their numerous generalizations represent one of successful mathematical models that 
^ \ is widely used in investigations of diverse physical phenomena and processes. 

For example, the following Euler's equation for an abstract model of an infinite-dimensional 
dissipative top 



j t A^^eB^ + [i>,A^]=ef, (1) 



where [ib, Aib] is Poisson's bracket, may be used to describe nonsteady-state flat-parallel flow of 
a viscous incompressible fluid in a channel with solid walls. The equation describing the fluid 
motion writes: 

d . . . . , dib dAib dib dAib , s 

--Aib + eAAib-^—^- + ^-—^ = e C osy, (2) 
at oy ox ox oy 

where ib(t,x,y) is a function of current, A = d 2 /dx 2 + d 2 /dy 2 is the Laplace operator. 
Under the following boundary conditions 

2-7T 

1) < x < — , < y < 2tt, a > 0, 

a 

2) ib (t,x + -^-,yj = ip(t,x,y), 



3) ib\ y=0 = — 



dib 
dy 



= ib\ = — 

y=0 ^ I 2/=2tt Qy 



y=2ir 
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(condition 3 indicates to the zero fluid flow rate) and the following correlation between the 
operators 

A S -A, B = AA, [*,„] = £3£_£«£, 

oy ox ox oy 

equation (J2J) is similar to (JJ). 

Such analogies allow one to conduct, for example, analysis of stability in problems of above 
type by classical methods of rigid body dynamics, and to suggest clear interpretation of results 
obtained [TUj . 

The present paper represents some results of qualitative analysis of the differential equations 
describing the motion of a rigid body in ideal incompressible fluid. If the following conditions 
are satisfied here, i.e. the fluid possesses a single-valued potential of rates and rests at infinity, 
then the body motion equations (6 ODEs) separate from the partial differential equations which 
describe the motion of fluid. In this case, the motion equations of the body coincide in their 
form with the corresponding Euler-Poisson equations and are called Kirchhoff's equations [HE]: 

. w dH dH dH 

M = MX 9M +7X 57' 7 = 7X 5M' 

where M = (Mi,M 2 ,M 3 ), 7 = (71,72,73) are vectors of "impulse moment" and "impulsive 
force", respectively. 

The total kinetic energy for the body and the fluid writes: 

2T = 2H = (AM, M) + 2(BM, 7) + (C7, 7). 

Here A, B, C are constant matrices. The latter are inertial characteristics of the body and the 
fluid. By a special choice of the origin and the direction of axes in the body it is possible to 
make the matrix A diagonal, and the matrix C symmetric, respectively. 

In problems of qualitative analysis of Euler's equations, it is possible to obtain the most 
complete results when the equations have many first integrals, for example, the equations are 
completely integrable. In this case, the phase space of system of the equations has a simple 
structure. 

We consider Kirchhoff's differential equations in Sokolov's case The equations in this 
case represent a completely integrable system. These write: 

Mi = M 2 M 3 + a(7 2 Mi + 7l M 2 ) + 2/3( 72 M 2 - 73M3) - 4 72 7 3 (2a 2 + 1 ) + 4a^7i7 3 , 
M 2 = 47i 73 (a 2 + 2/? 2 ) - /3(7 2 Mi + 71 M 2 ) - MiM 3 - 2a( 7 iMi - 73 M 3 ) - 4a^7 273 , 
M 3 = 4 7 i 72 (a 2 - 1 ) + /?( 73 Mi + 71M3) - a( 73 M 2 + 72 M 3 ) - 4a/3( 7l 2 - 7 2 2 ), 

71 = 72 (2M 3 + a 7 i) + /3( 72 2 - 7 f ) - 73 M 2 , 

72 = -7i(2M 3 + £72) - a(7? - ll) + 73M1, 

7 3 = 71 (M 2 + /?7 3 ) - 72 (Mi + a 73 ). (3) 

Here a, (3 are arbitrary constants. 

The system (J3J has the following 4 algebraic first integrals: 

2H = M 2 + Mf + 2M| + 2a( 73 Mi + 71 M 3 ) + 2/3( 73 M 2 + 72M3) + 4(071 - a-/ 2 ) 2 

- 4 73 2 (a 2 + (3 2 ) = 2h, 
V 1 = 7iMi + 7 2 M 2 + 7 3 M 3 = a, V 2 = 7 2 + 7I + 7 | = c 2 , 

^3 = {3(/?7i - arf2)(pMi - aM 2 ) + (2a 7 i + 2/3 72 + M 3 )((q 2 + /3 2 ) 73 + aM x + /3M 2 )} 2 
+ (M 3 - an ~ /?72) 2 {(/?Mi - aM 2 f + (a 2 + /3 2 )(2q 7 i + 2/3 72 + M 3 ) 2 } = c 3 . (4) 
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A lot of works were devoted to investigation of the Kirchhoff 's equations. A substantial part 
of these works is related to the problems of integrability. 

Another part of such works was devoted to the problems of investigation of stability for 
permanent motions (in particular, helical motions) of a rigid body in ideal fluid. The first 
results in this direction go back to Lyapunov [S]. This investigation was further developed, and 
some results can be found, for example, in [31 151 IS1IT3"]. 

The objective of the present work is to conduct qualitative analysis of solutions of equa- 
tions (J3J). We investigate a class of solutions of the equations, on which the elements of algebra 
of the problem's first integrals assume stationary values. Such solutions will be called stationary 
(see Appendix A for details). In particular, we found families of stationary solutions and families 
of invariant manifolds of steady motions (IMSMs) for the system of equations under scrutiny. 
We obtained conditions of stability and instability for several families of stationary solutions and 
families of IMSMs; parametric analysis of some of these conditions was conducted. Besides, some 
problems of bifurcations for both the families of stationary solutions and the families of IMSMs 
branching from these solutions (in particular, a trivial solution) were considered. Furthermore, 
the character of stability for branching manifolds was taken into account. 

The methods of investigation are based on classical Lyapunov's results (HI El, in particular, 
on his 2nd method. Methods of computer algebra implemented in CAS "Mathematica" were 
also used. A combination of these methods enabled to obtain results quite interesting from our 
viewpoint. 



2 Obtaining stationary solutions 

Let us consider the problem of finding stationary solutions and invariant manifolds of steady 
motions for the system ©. 

Analysis of equations (JHJ in terms of initial variables Mj, ji, i = 1,2,3, is rather bulky, and 
hence difficult. Therefore, in the papers devoted to the analysis of above equations, different 
linear transformations of variables are applied allowing one to reduce the equations and the inte- 
grals to a more compact form. In the present paper we use the following linear (not degenerate) 
transformation of the variables from pQ: 

Mt = st - ^dr 3 , M 2 = s 2 - ^/3r 3 , M 3 = s 3 + ^an + i/3r 2 , 

7i = r i; a = ^d, /3 = -ft, i = 1,2,3. (5) 

The latter enabled us to find out stationary solutions and IMSMs for the system (J2J) and to 
perform their analysis without going beyond standard algorithms. 

On account of the linear transformation of the variables ©, the equations of motion JHJ) for 
(3 = will take the form: 

h = (an + 2s 3 )r 2 - r 3 s 2 , si = (an + s 3 )s 2 - a 2 r 2 r 3 , 
r 2 = r 3 si - n(an + 2s 3 ), s 2 = (ar 3 - si)(an + s 3 ), 

r 3 = ns 2 - r 2 si, s 3 = -ar 2 s 3 , (6) 

and the corresponding first integrals write: 

2H = (si + s 2 2 + 2s 2 3 ) + 2ans 3 - a 2 r% = 2h, 

Vi = sin + s 2 r 2 + s 3 r 3 = a, V 2 = r\ + r\ + r\ = c 2 , 

2V 3 = (ansi + ar 2 s 2 + sis 3 ) 2 + s 3 (s| + (an + s 3 ) 2 ) = 2c 3 . (7) 



4 



V.D. Irtegov and T.N. Titorenko 



We shall now consider the problem of finding stationary solutions and IMSMs for the sys- 
tem (JSJ). We shall apply Routh-Lyapunov's method ^2] (see also Appendix B) for solving it. 
This will allow one to perform a substantial part of computational work with the use of computer 
algebra systems. 

In accordance with the method of Routh-Lyapunov, the functions K are constructed on 
the basis of the problem's first integrals. We shall construct here linear combinations of these 
integrals only (the combinations may be nonlinear ones): 

K = X H- - X 2 V 2 - X 3 V 3 , A; = const. (8) 

The integral K represents a family of first integrals, which is parametrized by the values Ao, Ai, 
A2 , A3 . We enter an excessive number of parameters into K that allows us to obtain "incomplete" 
combinations of integrals by equating some part of the parameters in K to zero. It is worth to 
note that different elements of algebra of first integrals will, generally speaking, correspond to 
various stationary solutions and manifolds. 

Next, we write down stationary conditions for K with respect to all the variables si, s 2 , s 3 , 

Ao«i - Airi - a 2 \ 3 rfsi - a 2 X 3 rir 2 s 2 - 2aA 3 riSis 3 - aX 3 r 2 s 2 s 3 - A 3 sis| = 0, 
X S2 - X\r 2 - a 2 X 3 rir 2 si - a 2 A 3 r|s 2 - aA 3 r 2 sis 3 - A 3 s 2 s| = 0, 
aA ri - Air 3 - aX 3 ris\ - aX 3 r 2 s 1 S2 + 2X s 3 - a 2 A 3 r 2 s 3 - A 3 sf s 3 - A 3 s|s 3 
- 3aX 3 r lS 2 3 - 2A 3 s| = 0, 

OOO Q 

o:Xqs 3 — 2A 2 n — A1S1 — a X 3 r\s 1 — a X 3 r2SiS2 — 0^3^53 — a A3ris 3 — CKA3S3 = 0, 
2A 2 r 2 + Ais 2 + a 2 A 3 riSis 2 + a 2 A 3 r 2 s 2 + aX 3 sis 2 s 3 = 0, 

(a 2 A + 2A 2 )r 3 + A1S3 = 0. (9) 

Solutions of system © define stationary solutions and IMSMs for the system (jHJ) (see Appen- 
dices A, B). In the general case these solutions may contain parameters Aj that appear in the 
family of integrals K (jHJ), and may hence represent a family of stationary solutions and IMSMs. 
Hence, to solve the problem formulated (obtaining stationary solutions and IMSMs for the sys- 
tem @ corresponding to the family of first integrals K) it is necessary to solve the system of 
6 algebraic equations containing four parameters Ao, Ai, A 2 , A3. The number of variables is 6. 
The system will have a closed form if any three first integrals (J7J are added to it. In particular, 
this means that the parameters Aj may be obtained as some functions of the constants of first 
integrals. 

All the equations of system Q are nonlinear, except for the last one that allows us to slightly 
simplify the problem. After removing the variable r 3 from the remaining equations @ with the 
use of the last one, we obtain the following system of 5 nonlinear algebraic equations containing 
5 variables: 

A si - Am - o 2 A 3 r 2 si - a 2 A 3 rir 2 s 2 - 2aA 3 riSis 3 - aA 3 r 2 s 2 s 3 - A 3 sis| = 0, 
A s 2 - Air 2 - a 2 A 3 rir 2 si - a 2 A 3 r|s 2 - aA 3 r 2 sis 3 - A 3 s 2 s| = 0, 
A 2 

Q AoH or — 7TT7-S3 - oiX 3 ris\ - aX 3 r 2 sis 2 + < 2.X Q s 3 - a 2 X 3 r\s 3 - X 3 s\s 3 

a z X + 2A 2 



ri, r 2 , r 3 : 
dK 

dK 

ds 2 
dK 

ds 3 

dK 
dr\ 
dK 

dr 2 
dK 
dr 3 
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- A 3 s 2 s 3 - 3aA 3 ris| - 2A3S3 = 0, 

OOO Q 

aAoS3 — 2A2?~i — A1S1 — a X3r±s 1 — a \3r2S1S2 — aX^s^s — a A3ris 3 — CKA3S3 = 0, 
2A 2 r 2 + Ais 2 + a 2 \ 3 risis 2 + a 2 A 3 r 2 S2 + aA 3 sis 2 S3 = 0. (10) 

The maximum degree of the equations belonging to the system is 3. 

To the end of obtaining solutions of system ()10|) . we apply the Grobner bases (GB) method [2] 
traditionally used in computer algebra for solving similar systems. Software implementation of 
the method can be found in many CAS. Application of only standard tools of CAS "Mathema- 
tica" for computing the GB allowed us to construct the Grobner basis for the system (|10|) under 
the following lexicographic ordering of the variables: r\ > r 2 > s 2 > s\ > S3. The timing for 
construction of the basis, as measured on a 1100 MHz Pentium with 256 MB RAM running under 
Windows XP, is 2.71 seconds. Below one can find a structure of Grobner's basis constructed 
under the indicated ranging of the variables for the given system of equations. 

S2S3/i(si, s 3 ) = 0, ((a 2 A + 2A 2 )si + aXis 3 )f 2 (si, s 3 ) = 0, S2/3O1, s 2 , s 3 ) = 0, 

S3/4(S1,S 3 )/ 5 (S 3 ) = 0, /6(^2,Sl,S 2 ,S 3 ) = 0, S 3 / 7 (si, S 2 , S 3 ) = 0, 

/sOi, s 2 ,s 3 ) = 0, /g(ri,si,s 3 ) = 0. (11) 

Here fi, i = 1, . . . , 9, are polynomials of the variables s±, s 2 , S3, r%, r 2 . The maximum degree of 
the polynomials is 7. These are too bulky, and so are omitted here. The system (|11|) is given in 
complete form in Appendix C. 

As is obvious from (|11|) . the basis constructed can easily be factorized that allows to de- 
compose it into several subsystems, which may be analyzed separately. Up to 12 subsystems 
were identified. For each of the subsystems we constructed a Grobner basis under lexicographic 
ordering of the variables. The latter enabled us to conduct some qualitative analysis of the set 
of solutions of each subsystem (with respect to the compatibility, finiteness or infiniteness of 
the set of the subsystems' solutions, etc.) and hence to obtain information about the whole set 
of system's (|11|) solutions (respectively, (fill]) and Q) and find out some groups of solutions. 

Finally, we conclude that the system © has an infinite set of solutions (the variable S3 is 
free). The following groups of solutions (besides the trivial solution) were found out: 4 families 
of IMSMs and 8 families of stationary solutions. Some of the solutions obtained can be found 
below. The solutions are given in the form representing the result of computing. 

1. The families of invariant manifolds of steady motions: 



si 



/2A2~(a 2 A + 2A 2 ) VA3S 2 - A + A X (2A 2 + a 2 A 3 sp 
a(a 2 X + 2A 2 )A 3 s 3 



v/2Aiii" ^/X 3 s 2 3 - A ((a 2 A + 2A 2 ) V / A 3 s^ - A + V^K 

S 2 = T / 9 . , M ! r l ~ 



a(a 2 Ao + 2A2)A3S3 ' a(a 2 Ao + 2A2)A3S3 

r 2 = t ;„ n . N , — , ^3 - 



a(a 2 A + 2A2)A 3 s 3 ' a 2 Ao + 2A2|' 



V2^(a 2 A + 2A 2 )v / A3s 2 ^Ao- A!(2A 2 + a 2 A 3 s 2 ) v / 2A 2 ^ 2 " 

s l = 7"^n 1 \ \ \ ) s 2 = T" 



a(a 2 A + 2A 2 )A 3 s 3 a(a 2 A + 2A 2 )A 3 s 3 ' 



v/A 3 si - A ((a 2 A + 2A 2 ) ^Ass^ - A - v^Ai) \/( x 34 - A )^ 2 

r i = — n \ u ) r 2 = ±- 



a(a 2 A + 2A 2 )A 3 s 3 ' a(a 2 A + 2A 2 )A 3 s 3 ' 

iTTSr)}- (12) 

a z X + 2X2 



1. 
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For brevity, we introduced the following denotations: 

Zl = a 4 X 3 + 4a 2 X 2 X 2 - 2XjX 2 + 4X X 2 2 - (a 4 A 2 , + 4A 2 . + a 2 {X\ + 4X X 2 ))X 3 s 2 3 



Z 2 



- 2 v / 2A 2 Ai(a 2 A + 2A 2 )^A 3 s 2 - A , 

a A Xl + 4a 2 X 2 X 2 - 2A 2 A 2 + 4A A^ - (a 4 A 2 , + 4A 2 , + a 2 {X\ + 4A A 2 ))A 3 s| 
+ 2y/2\^\ 1 (a 2 \ + 2A 2 ) J X 3 s 2 - A . 



2. The families of stationary solutions: 
< < n = 0, r 2 



V Ai\/— 2AqA 2 + 2AqA 2 



\/2a v / A 2 A3 



r 3 = 0, si = 0, 



S2 



VAi V /= 2A 2 " + 2A 2v % 
T ; . , «3 = 



a 



y— \A0A3 



'1 



»*3 



S3 



, ^-z 3 (2A 2 + ^ii)(a 2 A + 2A 2 + ^ 
±- " 7= , r 2 = 0, 



Ai^-z 3 (2A 2 + ^3") 



a ^3VA 3 



Ai^-2 3 (2A 2 + ^/i3") 
± — ^= , si 



az 3 y/X z 



s 2 =0, 



(a 2 A + 2A 2 ) V / -2 3 (2A 2 + ^) 



az 3 VA 



(13) 



Here 2 3 = a 4 A 2 , + 4A| + a 2 (A? + 4A A 2 ). 

Tools of computer algebra allow us to write rather easily the obtained solutions in terms of 
variables of any other "good" coordinate system, for example, in terms of the initial variables 
Mi, 7j. The solutions remain stationary also in terms of these variables. 

The families of stationary solutions (|13j) are given in terms of the variables Mi , M 2 , M 3 , 71 , 
72, 73: 



7i =0, 



72 



Mo 



V Ai-v/^2AoA 2 + 2AoA 2 
M 3 = 0}, 



73 = 0, Mi 



0. 



7i 



± 



v /Ai7^2A 2 " + 2 v / A^A 2 

3a v/AoXj 
V-Pi(2A 2 + v^I)(9a 2 A + 2A 2 + jfi) 



9a 2 piy/X~ 3 



72 = 0, 



73 



M 3 



Ai^-pi(2A 2 + v ^T) 
± ;== j Mi 



3apWX 3 



2Ai^-pi(2A 2 + v ^T) 



^ (18a 2 A + 4 A 2 - v ^T) A /-pi(2A 2 + ^ 



9api v A 3 




M 2 = 0, 



Here pi = 81a 4 A 2 , + 4A| + 9a 2 (A? + 4A A 2 ). 

Analysis of IMSMs (|12[) showed that, after transforming these expressions to the form, which 
does not contain the problem's variables under radicals, we obtain one family of IMSMs. The 
system of equalities (|12|) is a representation of the latter in different maps only. The equations, 
which define the family of IMSMs, write: 



a A 2 A 2 s| + 2ana' 3 AiA3SiS3 + X^sfs^ — 2aioA 2 A 3 S3 + 4anaAiA 2 A 3 sis 3 + 2agA 2 = 0, 



„2 \2„2„2 



On One Approach to Investigation of Mechanical Systems 



7 



4 4 X 4 2 4,^2 2\ \3 4 , a 2\2\2 4 a 2 2\ \2 2 

a xl a A 3 s 2 s 3 + 4a 11 aga A2A3S2S3 + 4a 9 A 2 A 3 s 3 — 4a 11 asa A2A3S2S3 



- 8a 4 AiA 3 sl + 4a§Al = 0, 
a ii^3 s 3 + 2a^ 1 aA|ris| + a^^A^r^sl — 2aj\3S^ — la^aXoX^riSz + agAo = 0, 
a 9^3 s 3 + 2ai 1 a 9 a 2 A|r 2 S3 - 2a 3 A|s3 + af^Xjr^sj - l^a^d 1 '\\r 2 s\ + a\\\s\ 



+ 2a\ l a s a 2 \Q\lr 2 sl - 2a 2 A A 3 S3 + OqAq = 0, 
ai±r 3 + Ais 3 = 0, 



(14) 



where a^, i = 1, . . . , 11, are polynomials of the parameters Ao, Ai, A2, A3. These are too cum- 
bersome, and so are omitted herein. 

From the geometric viewpoint, equations (|14|) - for each fixed collection of parameters of the 
family Ao, Ai, A2, A3 - describe the curves lying in M 6 at the intersection of the three 4th-order 
hypersurfaces, the 8th-order hypersurface and the hyperplane. 

The complete analysis of the family of IMSMs (|14|) and of motions on it is not given in the 
present paper. This analysis is rather nontrivial and may be a subject matter for another paper. 

3 Investigation of stability for stationary solutions 
and for IMSMs 

Let us consider the problem of stability for a series of families of stationary solutions and families 
of IMSMs of system ©. 

3.1 Investigation of invariant manifolds 

Consider equations © under the condition Ao = Ai = A3 = 0. It can be readily seen that thy 
have the following solution r\ = r 2 = r% = 0, which defines the 3-dimensional invariant manifold 
of steady motions for the equations @. The latter allow one to define the vector field on the 
IMSM (i.e. to reduce the initial system of equations to the IMSM): 

s'l = S2S3, s 2 = -S1S3, s' 3 = 0. (15) 

If we consider the initial values of the variables Sj, i = 1,2,3, as parameters, it is possible to 
assume that here we deal with the family of invariant manifolds. Hence, for S3 = this family 
adjuncts to the zero solution of the problem, i.e. it has at least one common point with the zero 
solution. 

When using the first integral V 2 = r\ + r\ + r|, which assumes the minimum value on the 
family of IMSMs, we easily conclude on stability of the elements of the family of IMSMs (|15|) . 
The vector field on the elements of the family of IMSMs has the two first integrals: 

W\ = s\ + s| = mi, W 2 = s 3 = m 2 . 

Consequently, in this case we may speak of investigation of the 2nd-level stationary solutions 
(see Appendix A). Finding such solutions and their analysis are trivial in the computational 
aspect, and we will not concentrate on them here. Let us consider a more complex case. 

3.2 Investigation of 2nd-level stationary solutions 



By constructing the Grobner bases with respect to the problem's variables and to one or several 
parameters Aq, Ai, A2, A3 for the subsystems of system ((TTj) we can also obtain solutions under 
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some conditions imposed on above parameters. For example, we constructed the Grobner basis 
with respect to the variables n, r 2 , s\ and the parameter A2 for the subsystem 

S2S3 = 0, ((a 2 A + 2A 2 )si + aAis 3 ) / 2 (si,s 3 ) = 0, s 2 = 0, s 3 = 0, 
/6(^2,si,s 2 ,S3) = 0, /80i,S2,s 3 ) = 0, /g(n,si,s 3 ) = (16) 

of system IJ11J1 . It enabled us to obtain solutions of the system under the following conditions 
imposed on the parameter A 2 : 

2a 2 A 2 + A 2 + A 1 y4a 2 A ! + A 2 h7 , 
X2 = 4An" ' (17) 



2a 2 A 2 + A 2 -A 1 V4a 2 A 2 + A 2 
A2 " 4A^ • (18) 

The solutions obtained can be found in Appendix D. The solutions are given right in the form 
representing the result of computing. 

Likewise in case of Q12|). after transforming the expressions of solutions obtained to the form, 
which does not contain the problem's variables under radicals, we found out that these solutions 
represent the families of IMSMs for the system © , these solutions being written in terms of the 
maps for these families. Finally, we have 2 families of IMSMs. 

The first family of IMSMs, which corresponds to A 2 (|17|). can be written as: 



2a 4 X 2 r 2 + (^2a 2 X 2 + A x f Ai - ^4a 2 A 2 + A 2 ^ 
\a 2 X 2 + Ai (Ai - ^4a 2 A 2 + A 2 ^ 



si 



Ao 

A3' 



si + 4 = s 3 = 0, r 3 = 0, 2a 2 A r 2 + ^Ai - ^4a 2 A 2 + Xj^j s 2 = 0, (19) 

and the second family of IMSMs, which corresponds to A 2 (|18|). can be written as: 
2a 4 X 2 r 2 + (2a 2 Xl + Ai [Xi + ^4a 2 A 2 + A 2 ^ 
'2a 2 A 2 + X X U + ^4a 2 A 2 + A 2 )) ^, 



4 



4 + 4 = ^, s 3 = 0, r 3 = 0, 2a 2 A r 2 + ^Ai + ^4a 2 A 2 + A?) s 2 = 0. (20) 



A 



3 

From the geometric viewpoint, each of these families of IMSMs - for each fixed collection of 
parameters of the family Ao, Ai, A 3 - describes the curves lying in R 6 at the intersection of the 
three hyperplanes, one elliptic "cylinder" and one circular "cylinder". 

The vector field on elements of the family of IMSMs (|19[) is given by the differential equation: 

■ = 2aAo(A - A3S 2 ,) 

(A 1 + V4a 2 A 2 + A 2 )A 3 ' 1 j 

which is derived from equations © after removing s±, s 3 , r±, r 2 , r 3 from them with use of 
expressions (fTU|) . 

The vector field on elements of the family of IMSMs (|2U|) is given by the differential equation: 

2aA (Ao - A3S 2 .) 
(Ai - X /5?A 2 + A 2 )A 3 
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which can also be derived from equations (jHJ) after removing si, S3, r% from them with 

the use of (j2U|). 

Consider the problem of finding 2nd-level stationary solutions for the system ()21|) and inves- 
tigation of their stability on elements of the family of IMSMs (jl9|) . 
As obvious from 1)21(1 , solutions of the form 

^2 = 1 — 5 ^2 = — 1 — (^*^) 

V A3 V A3 



are its stationary solutions. Now we investigate their stability by Lyapunov's method [Sj. Let 
us consider the first of these solutions (J2*H|) . 

According to the above method, we consider the solution s\ = — \/Ao/V^3 m the capacity 
of undisturbed one. Next, we introduce the deviations z = s 2 — s® of disturbed motion from 
undisturbed one, and consider the function V = \z 2 in the capacity of the Lyapunov function. 

The equation of disturbed motion writes: 

. _ 2a\ z 2 Aa\l /2 z ^ 



Ai + v/4a 2 Ag + A| (Ai + V4a 2 Ag + A=f)A 3 

The derivative of the function V due to the differential equation (|24() up to the 2nd order terms 
writes: 

4aA 3/2 z 2 
V = + V 

(A 1 + v /4a^ + A 2 )A 3 

where V n is the n-order terms (n > 2). 

The function is positive-definite. Prom the expression for its derivative it is obvious that 
it will be negative-definite when the following conditions imposed on the parameters A, hold: 

a > V A > V (Ai < V A 3 > A Ai > V A 3 > 0). (25) 

According to the Lyapunov theorem [§] on stability of undisturbed motion, satisfaction of these 
conditions means that the solution investigated is asymptotically stable on elements of the family 
of IMSMs P]). 

According to the above theorem, the solution under scrutiny is unstable on elements of the 
family of IMSMs (|19|) when the following conditions hold: 

a < V A > V (Ai < V A 3 > A Ai > V A 3 > 0). (26) 

Substitution of the stationary solution s 2 = — V^o/V^? into allows one to obtain the 
solution corresponding to it in the whole space of variables s\, S2, S3, ri, r 2 , ^3: 



< n = 0, r 2 = v _ " -, r 3 = 0, si = 0, s 2 = ■==, s 3 = 0} . (27) 

[ 2a^VAoA3 VA3 J 

Analysis of (|27|) showed that this solution is unstable in the sense of Lyapunov [H] under the 
following conditions imposed on the parameters A^: 

a<OV(Ao<OVAi>OVA 3 <OAAo>OVAi>OVA 3 >0) 

A a >0 V (A < V Ai >0 V A 3 < A A > V Ai > V A 3 >0). (28) 

There are roots having positive real part among the roots of the characteristic equation con- 
structed for the equations © linearized in the neighbourhood of the solution (|27|) when condi- 
tions (PHI) hold. 
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Comparison of the conditions (|25|) and (|28|) shows that under the same conditions imposed 
on the parameters Aj: 

a > V A > V Ai > V A 3 > 0, 

the first of the solutions (|23j) asymptotically stable on the IMSM (|19j) corresponds to the solution 
which is unstable in the whole space of the problem's variables. 

Similarly, comparison of the conditions Q26|) and 1)28(1 allows us to conclude: when the following 
conditions 

a < V A > V Ai > V A 3 > 

hold, the solution unstable on the IMSMs (|19|) corresponds to the solution unstable in the whole 
space of the problem's variables. 

We also investigated stability of the 2nd solution of (|23|) and of stationary solutions of sys- 
tem (|22() . The results appeared to be similar to those given above. 

As far as behaviour of the IMSMs themselves (|19|) . (|20j) is concerned, note the following. The 
manifolds intersect and coincide under the condition Ao = Ai = 0. As a result, we obtain the 
invariant manifold defined by the equations: 

81 = 0, s 2 = 0, s 3 = 0, r 3 = 0. (29) 

The vector field on IMSM (|29|) is given by the differential equations: 

f\ = arir2, r2 = —o?r\. 

4 On stability of the zero solution and of IMSMs adjunct to it 

Let us consider a class of special stationary solutions. These solutions possess the following 
properties: several of the problem's first integrals assume stationary values on the solutions 
and, as a rule, there is a bifurcation of stationary solutions (invariant manifolds) of various 
dimensions in their neighbourhood. 

For example, in completely integrable cases of the problem of rigid body's motion having 
one fixed point, when all the first integrals are quadratic, the manifolds of dimension 3 (it is 
half of the number of variables used in describing the problem) are typical IMSMs branching 
from special permanent rotations. Furthermore, in many cases there is the following relationship 
between the stability of special permanent rotations and the stability of IMSMs branching from 
them: branching of stable IMSMs from special permanent rotations represents the necessary 
and sufficient stability condition for them. If algebraic first integrals of the problem under 
consideration are not only quadratic then the relationship between the property of stability of 
special stationary solutions and the property of stability of IMSMs branching from them is more 
complex. 

As far as Kirchhoff 's equations are concerned that is obvious from their analysis in Sokolov's 
case, the IMSMs of both even and odd dimensions can branch from special stationary solutions 
(helical motions). Now, we consider a particular example of such a bifurcation of stationary 
solutions. 

Consider the problem of stability of both the zero solution and the IMSMs adjunct to it, i.e. 
the stability of invariant manifolds which have at least one common point with the zero solution. 
Equality of the system's © Jacobian, which is computed for zero values of the variables Sj, ri 
i = 1, . . . , 3, to zero is the condition of existence of such IMSMs. 

The Jacobian of © for s\ = 0, S2 = 0, S3 = 0, r\ = 0, r2 = 0, r^, = writes: 

J = (Xf + 2A A 2 ) {a 2 \l + \\ + 2A A 2 ) {a 2 X 2 + Xl + 4A A 2 ). 
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The following possibilities for satisfaction of the equality J = were considered (when the zero 
solution is special): 

A 2 

1)A = 0, Ai=0; 2)Ai = 0, A 2 = 0; 3) A 2 = -— L ; 

4)A 2 = -^i; 5)X2 = -^±A. (30) 

Obtaining solutions of stationary equations © for the indicated values of A 2 (|3U|) with the aid 
of Grobner's bases technique allowed us to find out a sufficiently large number of families of 
IMSMs adjunct to the zero solution. Some of these solutions are adduced below: 

Ai s 3 A \ a 2 A§ + A? 

si = r-s 3 , s 2 = 0, n = , r 2 = 0, r 3 = — s 3 } for A 2 = ^ ; (31) 

aAo a Ai J 2Ao 

{s 3 = 0, n = 0, r 2 = 0, r 3 = 0} for A = 0, Ai = 0; 

{si = 0, s 2 = 0, s 3 = 0, n = 0, r 3 = 0} for Ai =0, A 2 = 0; 

A 1 A 2 

si = 0, s 2 = — r 2 , s 3 = 0, n = 0, r 2 = ^ for A 2 = -777-, A 3 = 0; 

2aAoAi 2aAn 2aAoAi 

51 = -^aFa! S3 ' 52 = °' ri = -^Rf S3 ' r2 = °' r3 = "^aFa! S3 

n 2\ 2 1 \2 

for A 2 = - a 7 + Al , A 3 = 0; (32) 
4Ao 

jsi = 0, s 2 = 0, s 3 = ~rx, r 2 = 0, r 3 = o| 

^2x2 , \2 

for A 2 = - Q A ; + \ A 1= 0, A 3 = 0. 

4Ao 

These solutions represent the simplest of the invariant manifolds obtained: the zero vector field 
is defined on all the IMSMs (|31j) . (|32j) . The latter attributes some specificity to the problem of 
stability analysis for such manifolds. 

We use the method of Lyapunov functions, in particular, the Routh-Lyapunov method, to 
investigate stability of the zero solution as well as the families of IMSMs (|31j). (|3*2*|) adjunct to it. 

4.1 Investigation of stability of the zero solution 

The procedure of obtaining sufficient conditions of stability for the stationary solutions by 
Routh-Lyapunov 's method is practically reduced to the verification (in the simplest case) of 
signdefiniteness of the 2nd variation of the integral K (jHJ) in the neighbourhood of the statio- 
nary solution under scrutiny. The sufficient conditions can be made "softer" , if signdefiniteness 
of 8 2 K is considered on the manifold defined by the first variations of each of m — 1 integrals 
(where m is the number of vanishing integrals in K). 

Obtaining sufficient stability conditions for the zero solution by this technique is rather 
trivial in the computational aspect. Note only that it is stable in the sense of Lyapunov when 
the following restrictions imposed on the problem's parameters are satisfied: 

a>0VA o >0vfA 1 >QVA 2 <- a2A ^ +A ' AA 1 <0VA 2 <- a2A j, + A n A 
\ ^Ao 2Ao / 

a <0VA >0vfA 1 >0VA 2 <- a2A ° x +A ' AA 1< 0VA 2 <- q2A ° x + X " ) . (33) 
V ^Ao 2Ao / 

It is possible to slightly "soften" the stability conditions obtained, considering that the zero 
solution is special. To this end, we solve the problem of choosing the "best" first integral, i.e. 
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the one which gives the "most soft" sufficient stability conditions. The solution of this problem 
can be obtained, for example, with the use of parametric analysis of stability conditions, i.e. 
their minimization with respect to one or several parameters. 

Consider one of the conditions imposed on the parameters (I33|) : 

A2< _a2Af + A? or A = 2A2Ao + a 2 A 2_ A 2 <a (34) 

zAo 

Next, we find a stationary value A with respect to Ao: 

— = 2A 2 + 2a 2 A = 0. 
d\ 

The latter holds for Ao = — A 2 /a 2 . Its substitution into (|34(1 gives the following condition 
imposed on the parameters: 

\2 

-Va?<o. 

Hence we have stability of the equilibrium position under and Ai 7^ or a 7^ and A2 7^ 0. 

4.2 Investigation of stability of IMSMs adjunct to the zero solution 

Let us investigate stability of the family of IMSMs (|31|) . The vector field on the elements of the 
family of IMSMs (j31|l is given by the differential equation 

s 3 = 0, (35) 

derived from equations © after removing s±, s 2 , rx, r 2 , ^3 from them with the aid of expres- 
sions (|3"Tj) . 

Using 1)35(1 . we may conclude that the elements of the family of IMSMs obtained represent 
some curves in R 6 , over each point of which the one-dimensional family of solutions (S3 = S3 = 
const) for the equation (|35|) is defined. Such IMSMs with the "bundle" defined on them will 
be called "framed invariant manifolds". Each point in the framed IMSM corresponds to some 
helical motion of a rigid body. 

To the end of obtaining sufficient stability conditions for framed IMSMs we use the standard 
Lyapunov's technique. 

The second variation of K in the neighbourhood of some helical motion S3, which lies on the 
chosen IMSM represented in terms of deviations 

. S3 A s 3 Ais 3 

zi=n-\ , z 2 = r 2 , z 3 = r 3 — , z 4 = si H — , 

a Ai aAo 

Z5 = S 2 , Zq = S3 - S§ 

writes: 

(a 2 Ag + Af)(Ao-A 3 sf) q^g + \\ \\ A 2 

5 K ~ V? Zl + 2A " 2 + VT 3 " 114 + I 4 

- Aiz 2 z 5 + -(A - \ 3 sf)z%. 

And the respective variations of the first integrals H,V\,V2 are: 
a 2 An Ai Af - a 4 A 4 



- -0 A l a, - U An \ n 

6H = [azi — z 3 — z 4 H ? . 2 . 2 U Z6 U3 = °> 

aA a 2 AgAf 
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SVi = —zi + z 3 - -z 4 + y 2 .° . 1J z 6 s° = 

\ oaq a a^AoAi / 

SV 2 = 2 (-U + ^z 3 + ^f^M s° = 0. 
\ a Ai a z \f ) 



(36) 



Conditions of signdefiniteness for 5 2 K are sufficient stability conditions for the elements of the 
family of IMSMs ()31[) . After trivial transformations, these conditions may be written as follows: 

A - X 3 sf > 0, a 2 A[] - (a 2 A 2 , + A 2 ) \ 3 sf > 0, 

A? A 3 sf + a 2 A 2 Ai (\ 3 sf - A ) / 0, A > 0, A x ^ 0. 

Now we extend the problem and investigate stability of a body's helical motions corresponding 
to the elements of the framed IMSM. 

Assuming that s!] ^ 0, we eliminate the variables z\, Z4 with the use of equations (|36j) 
(because there are only two linearly independent ones) from 5 2 K. As a result, the following 
quadratic form yields: 

2 ~ a 2 A 2 + A 2 2 A| + a 2 Ao(a 2 A 2 + A 2 )(Ao-A 3g f) 2 1 2 2 
2A^ 22 2A A 2 3 ~ 225 2 ^ ~ 353 ' H 

(a 2 X 2 + X 2 )(a 2 Xl-(a 2 X 2 + X 2 )X 3 sf) 
+ X«X\ Z3Z6 

(a 2 X 2 + X 2 ) 2 (a 2 Xl-(a 2 X 2 + X 2 )X 3 sf) 2 
+ 2a 2 A 2 Af 26 • 

Conditions of signdefiniteness for 5 2 K are sufficient stability ones for the helical motions, 
which belong to IMSMs under scrutiny. When representing them in the form of the Sylvester 
conditions, we have: 



i) a2x \ +x Ko, 

' 2A 



(a 2 A 2 + A 2 )(AK« 2 Ao(a 2 A 2 + A 2 )(Ao-A 3g n) > n 
' 4A 2 A? > U ' 



3 ) ^r^(« 2 A|]-(a 2 A 2 + A 2 )A 3S f)(Af + a 2 A (a 2 A 2 + A 2 )(A -A 3S f)) >0 



8X X 1 



(a 2 A 2 + A 2 ) 3 (A 2 A 3S f + a 2 A 2 (A 3g f - A )) : 
16a 2 A^Af 



4) K - ' 11 y 1 6 3 ' 2^ >0 (37) 



A standard software package "Algebra InequalitySolve" of CAS "Mathematica" was used for the 
purpose of verification of compatibility for this system of inequalities. Its application to (|H7|) 
showed that the inequalities are compatible when: 



a < V A > 



\3/2 x 3/2 



VlA 1 <0V(A 3 <0AA 3 >0V ~\° = < s 3 < 



V(« 2 A 2 + A 2 )A 3 V(« 2 A^ + A 2 )A 3/ 

,3/2 x 3/2 



AAi > V I A 3 < A A 3 > V '° „ = < s\ < „° , == 

V V(« 2 A^ + A 2 )A 3 3 V(« 2 A^ + A 2 )A 3 , 
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Aa > V A > 

V( X 1 <Qy(x 3 <QyX 3 >0A—= ^ = < s °< - ^ =) 
V V V(« 2 A^ + A^A 3 V(« 2 A^ + ADA 3 y 

AA!>0vfA3<0VA 3 >0A ; ?$1 = < s° 3 < - a ^ —=\). 

V VV*8 + A?)A 3 3 v/(a^ + A?)A 3 J J 

Hence, stable in the sense of Lyapunov are only those body's helical motions for which the 
parameter satisfies the latter conditions. 

Comparison of stability conditions for the zero solution with the stability conditions for 
IMSMs 1)3 ljl gives evidence that satisfaction of the former conditions implies satisfaction of the 
latter. Consequently, the stable IMSMs adjunct to the stable zero solution. 

We also investigated stability of IMSMs 1)32)1 . We proved that some of these IMSMs are 
unstable with respect to the first approximation. Investigation of stability for other solutions 
necessitates involving of higher-order terms in the integral K expansion. But this problem was 
not considered here. 



5 Conclusion 

The paper represented some results of qualitative analysis of the differential equations describing 
the motion of a rigid body in ideal fluid. The Kirchhoff 's differential equations were considered 
in Sokolov's case. In this case, these have 4 algebraic first integrals (three quadratic ones and 
one 4th degree integral) and represent a completely integrable system. 

Routh-Lyapunov's method was used to analyze the set of solutions of the equations. The 
method proposes a technique of finding both the stationary solutions and the invariant manifolds 
of motion equations, when the equations have a sufficiently large number of first integrals. The 
stationary solutions and invariant manifolds, which were obtained by this method, may be 
investigated for stability by the 2nd Lyapunov's method. Furthermore, the corresponding first 
integrals can be used here as Lyapunov functions. 

The paper considers a rather typical case when the problem of finding stationary solutions is 
reduced to solving a nonlinear system of equations. Some interesting cases of stability investiga- 
tion for conservative systems are also given. For example, the property of asymptotic stability 
of the equilibrium state for a vector field on the 1-dimensional invariant manifold of a conserva- 
tive system was used for proving instability of the stationary solution (which corresponds to this 
equilibrium state) within the whole system's phase space. Some examples of branching invariant 
manifolds of various dimensions and investigation of their stability were considered. 

The results of this work were obtained with the use of computer algebra tools. Such investi- 
gations cannot probably be conducted within an acceptable time without them. To ground this 
statement, we gave a Grobner basis in Appendix C. This basis was constructed and used for 
finding solutions of a nonlinear system of algebraic equations arising in computations. 

The results of qualitative analysis of Kirchhoff's differential equations represented in the 
paper give evidence that the technique of investigation of mechanical systems, which is based 
on a combination of classical methods of rigid body dynamics and computer algebra methods, 
is rather efficient and may be used for investigations of above type problems. 



A Appendix 



The following concepts were used in the paper. 
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Definition 1. The solutions of differential equations, on which the first integral (an element 
of the algebra of problem's first integrals) assumes a stationary value, are called the stationary 
solutions. 

This means the following. The stationary solution satisfies the equations which are the 
result of equating all the partial derivatives of the first integral with respect to the problem's 
variables to zero (stationary conditions). Such solutions simultaneously are the solutions of 
initial differential equations. This follows from the Lyapunov's theorem j^]. 

Definition 2. Manifolds, whose equations satisfy the stationary conditions for some first inte- 
gral, are called the invariant manifolds of steady motions (IMSMs). 

The proof of invariance of such manifolds with respect to initial differential equations can be 
found in Appendix B. Note, that manifolds are understood as the sets of dimension larger than 
zero. 

If some differential equations have an IMSM then by these equations it is possible to define 
a vector field on this IMSM. This procedure is called a reduction of the initial system. The 
vector field, in turn, can have first integrals which may be used for finding stationary solutions 
and the 2nd-level IMSM (on the given IMSM). 

Stationary solutions and manifolds are suitable in stability investigations by Lyapunov's 
second method, because the expansions of the corresponding first integrals in the neighbourhood 
of these solutions and manifolds in Taylor series do not contain linear terms. 



B Appendix 

More exactly, we use the following theorem: 

Theorem 1. If partial derivatives of the first integral V(x,t) of the system Xi = Xi(x,t) (i = 
1, . . . ,n) with respect to the problem's phase variables have the form 

dV ^ k k 

— = ^ a il( X > ^Vifo *) + O^O^ *)Pp( x > *H » i = 1, ... ,71, 

X% 1=1 1=1 p=l 

and the rank of the matrix \\au(x, t)\\ is "k" on the manifold ipi(x, t) = 0, I = 1, . . . , k, then the 
manifold <pi(x,t) = 0, I = 1, . . . , k, is invariant for the initial system of differential equations. 

Proof. Let the system of differential equations 

X{ — A^(x, £), i — 1,2,..., 7i, 
have the first integral V(x, t) and 

dV fc ^ k r fc 

= ^2 a u{x,t)tpi(x,t) + ^2^2a Up (x,t)tpi(x,t)tp p (x,t) H , i = l,...,n. 

% i=i i=i p=i 

Since V(x,t) is the first integral, we have 

~dt ~ ^ dx~ j + ~dt~ 

When differentiating the latter identity with respect to Xi, we obtain the system of equations 

A d fdv\ x 8 8V ^ dVdXj _ Q f _ l 

^— ' cIt.; \ f)nr j I 3 Bt.: Fit QXj QXi 

3=1 



J dxi \ dxj J dxi dt 
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Having changed the order of differentiation in the above system, we have 




dV dXj 

dxj dxi 



i = 1, . . . , n. 



Now substitute the expression for dV/dxi into the latter formulae and perform differentiation 
with respect to x^. After trivial transformations we have: 



where F;(0) = 0. 

When ip p = 0, p = 1, . . . , k, the latter system transforms into the system of n linear homoge- 
neous equations: 



Under the condition that the rank of the matrix ||aj/(a;,f)|| is k on the manifold ip p = 0, p = 
1, . . . , k, this system has only the following trivial solution: 



The latter proves that the manifold (p p = 0, p = 1, . . . , k is invariant for the initial system of 



C Appendix 

The Grobner basis constructed for the system pui) : 

s 2 S3{2A 2 (a 4 A[j + W\l\ 2 + 2A 2 (A 2 + 2A A 2 )) + 4aAiA 2 (a 2 A + 2A 2 )A 3 sis 3 
- 2A 2 (a 4 A 2 + 4A 2 - 2a 2 {X\ - 2A A 2 ))A 3 s 2 + a 2 (a 2 X + 2X 2 ) 2 X 2 s 2 s 2 
+ 2a 3 Ai(a 2 A + 2A 2 )A 2 ; s 1 s 3 i + cZA^A 2 ^} = 0, 

((a 2 A + 2A 2 )si + aAis 3 ){2A 2 (a 4 Aij + 4a 2 X 2 X 2 + 2A 2 (A 2 + 2A A 2 )) 

+ 4aAiA 2 (a 2 A + 2A 2 )A 3 sis 3 - 2A 2 (a 4 A 2 1 + 4A 2 , - 2a 2 (A 2 - 2A A 2 ))A 3 s 3 ! 
+ a 2 (a 2 A + 2A 2 ) 2 Als 2 sl + 2a 3 Ai(a 2 A + 2A 2 )A^is| + a 4 A 2 Als 4 } = 0, 

s 2 {2A 2 (a 2 A + 2A 2 ) 3 (A 2 + 2A A 2 ) + 4a 2 A A 2 (a 2 A + 2A 2 )(a 2 A - aAi + 2A 2 ) 

x (a(aX + Ai) + 2A 2 )A 3 s 2 - a 4 A (a 2 A + 2A 2 ) 3 Als 4 + 4a 2 A A 2 (a 2 A + 2X 2 fX 3 s 2 2 

+ a 4 A (a 2 A + 2A 2 ) 3 A 2 ,s 2 1 + 4aAiA 2 (2a 6 A 3 1 + 16a 2 AoA 2 + 8A 3 , 

+ a 4 A (10A A 2 - A 2 ))A 3 sis 3 - 4a 5 A Ai(a 2 A + 2A 2 ) 2 A!s?s 3 - 2A 2 (a 2 A + 2A 2 ) 

x (8Al - 4a 2 A 2 (A 2 - 2A A 2 ) + a 4 A (2A A 2 - A 2 ))A 3 sl + a 2 (a 2 A + 2A 2 ) 

x (8a 2 A x X 2 2 + 8A 3 , + a 4 A (-5A 2 + 2A A 2 ))A|s 2 s| + 2a 3 Ai(8a 2 A A 2 , + 8A 3 , 

x +a 4 A (2A A 2 - A 2 ))A^is| + 2a 4 A 2 A 2 (a 2 A + 2A 2 )A|s 4 } = 0, 

s 3 {2X 2 (a 4 X 3 + 4a 2 A§A 2 + 2A 2 (A 2 + 2A A 2 )) + 4aAiA 2 (a 2 A + 2A 2 )A 3 sis 3 






differential equations. 
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- 2A 2 (a 4 A 2 1 + 4A 2 , - 2a 2 {Xj - 2A A 2 ))A 3 s| + a 2 (a 2 A + 2A 2 ) 2 A|s?s| 

+ 2a 3 Ai(a 2 A + 2A 2 )A| Sl si + a 4 A 2 A|s|}{(a 2 A + 2A 2 ) 4 (a 2 A 2 + A 2 + 4A A 2 ) 

- 4A 2 (a 2 A + 2A 2 ) 2 (a 4 A 2 + 4A 2 + a 2 (A 2 + 4A A 2 ))A 3 s 2 - a 2 (a 4 A 2 + 4A 2 
+ a 2 (A 2 + 4A A 2 )) 2 A 2 s 4 } =0, 

2AiA 2 (a 2 A + 2A 2 ) 2 r 2 - 2A A 2 (a 2 A + 2A 2 ) 2 s 2 - A (a 3 A + 2aA 2 ) 2 A 3 s 2 s 2 

- A (a 3 A + 2aA 2 ) 2 A 3 s 2 i - aAi(a 2 A + 2A 2 ) 2 A 3 sis 2 s 3 + 2A 2 (a 2 A - aAi + 2A 2 ) 
x (a(aX + Ai) + 2A 2 )A 3 s 2 s| - a 2 (a 2 A + 2A 2 ) 2 A|s 2 s 2 s| 

- 2a 3 Ai(a 2 A + 2A 2 )A 2 .s :L s 2 s 3 ! - a 4 A 2 A|s 2 s| = 0, 

s 3 { - 2A 2 (a 2 A + 2A 2 ) 8 (a 2 A 2 1 + A 2 + 4A A 2 ) + (a 2 A + 2A 2 ) 2 (a 16 A 8 1 + 256A 8 , 
+ 16a 12 A^A 2 (A 2 + 7A A 2 ) + a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A^A 2 + 16A A 2 ) 
+ 32a 4 A|(A 4 + 16A A 2 A 2 + 56AqA 2 ) + 4a 10 A 3 ,A 2 (A 4 + 25A A 2 A 2 + 112X 2 X 2 2 ) 
+ 8a 8 AgA^(3A 4 + 40A A 2 A 2 + 140A^A|) + 4a 6 A|(A? + 12A A 4 A 2 + 140A^A 2 A| 
+ 448A 3 A 3 ))A 3 s 2 + (a 2 A + 2A 2 ) 2 (a 16 A 8 + 256A 8 + 16a 12 A^A 2 (A 2 + 7A A 2 ) 
+ a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A^A 2 + 16A A 2 ) + 32a 4 A^(A 4 + 16A A 2 A 2 + 56X 2 X 2 2 ) 
+ 4a 10 A 3 1 A 2 (A 4 + 25A A 2 A 2 + 112AgA|) + 8a 8 A^A^(3A 4 + 40A A 2 A 2 + 140A^A|) 
+ 4a 6 (A? + 12A A 4 A 2 + 140AgA 2 A^ + 448A 3 1 A|))A 3 s 2 , + 2aAi(a 2 A + 2A 2 ) 3 
x (a 12 X 6 - 96a 2 A A^ - 64X 6 2 - 16a 4 A A^(A A 2 - A 2 ) + 2a 6 A A 2 (A 2 + 4A A 2 ) 
x (A 2 + 6A A 2 ) + 4a 8 AjjA 2 (2A? + 9A A 2 ) + a 10 A 4 (A 2 + 10A A 2 ))A 3 sis 3 
+ (a 18 A^A 2 + 2048A^0 - 256a 2 A^A 2 - 32A A 2 ) + a 16 A^(A 4 + 12A A 2 A 2 + 8X 2 Xl) 
+ 8a 14 A^A 2 (A 4 + 7A A 2 A 2 + l6X 2 Xl) + 64a 4 A 2 1 (224A 2 A 2 , - 3A 4 - 12A A 2 A 2 ) 
+ 4a 12 A 4 A 2 ,(3A 4 + 28A A 2 A 2 + 224A 2 1 A 2 ,) - 32a 6 Af(A? + 12A A 4 A 2 + 28A 2 ,A 2 A 2 . 

- 448AgA|) + 8a 10 X 2 X 2 2 {U8XlX 3 2 - A? - 8A A 4 A 2 ) - 4a 8 A^(A? + 8A A?A 2 + 68X 2 XjX 2 2 
+ 112A;]A 2 A 3 - 2240A 4 A|))A 3 s|, - a 2 (a 2 A + 2A 2 ) 2 (a 4 A 2 1 + 4A 2 + a 2 (A 2 + 4A A 2 )) 

x (a 10 A[] + 6AX 5 2 + a 8 Aj](A 2 + 12A A 2 ) + 8a 2 A^(A 2 + 18A A 2 ) + 2a 6 A^A 2 (3A 2 

+ 28A A 2 ) + 2a 4 A 2 (A 4 + 6A A 2 A 2 + 64A 2 ,A 2 ,))A 2 ; s 2 s| - 2a 3 Ai(a 2 A + 2A 2 )(a 4 A 2 1 + 4A| 

+ a 2 (A 2 + 4A A 2 ))(a 10 A^ + 96X 5 2 + 8a 6 AgA 2 (A 2 + 9A A 2 ) + 16a 2 A^A 2 + 13A A 2 ) 

+ a 8 A 3 (A 2 + 14A A 2 ) + 2a 4 A 2 (A 4 + 10A A 2 A 2 + 88A 2 A 2 ))A 2 s lS 3 + a 2 (a 4 A 2 + 4A 2 

+ a 2 (A 2 + 4A A 2 ))(128A 2 7 - 96a 2 A|(A? - 4A A 2 ) + a 12 A^(2A A 2 - A 2 ) 

+ a 10 Xl{24X 2 X 2 2 - Xj - 14A A 2 A 2 ) + 2a 8 A 2 1 A 2 (60A 2 l A 2 , - 5Xj - 36A A 2 A 2 ) 

+ 8a 4 A|(60A^A| - 3A 4 - 26A A 2 A 2 ) - 2a 6 A 2 (A? + 14A A 4 A 2 + 88X 2 XfX 2 2 

- 160A 3 A 3 ))A|s 4 - (a 3 Ac + 2aA 2 ) 4 (a 4 A 2 + 4A 2 + a 2 (A? + 4A A 2 )) 2 A 3 s 2 S 4 

- 2a 5 Ai(a 2 A + 2A 2 ) 3 (a 4 A 2 1 + 4A 2 , + a 2 {Xj + 4A A 2 )) 2 A 3 ! sis| - a 6 A 2 (a 2 A + 2A 2 ) 2 
x (a 4 A 2 + 4A 2 + a 2 (A 2 + 4A A 2 )) 2 A 3 ^} = 0, 

-2A 2 A 2 (a 2 A + 2A 2 )(a 16 A 8 + 256A 8 + 16a 12 A^A 2 (A 2 + 7A A 2 ) + a 14 A^(A 2 + 16A A 2 ) 
+ 64a 2 A^(3A 2 + 16A A 2 ) + 32a 4 A^(A 4 + 16A A 2 A 2 + 56AgA|) 
+ 4a 10 A 3 1 A 2 (A 4 + 25A A 2 A 2 + 112AgA|) + 8a 8 A^A^(3A 4 + 40A A 2 A 2 + 140AgA|) 
+ 4a 6 A^(Af + 12A A 4 A 2 + 140A§A?A| + 448A 3 ,A 3 ,))si - (a 2 A + 2A 2 ) 3 (a 16 A 8 1 + 256A 8 , 
+ 16a 12 A^A 2 (A 2 + 7A A 2 ) + a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A^A 2 + 16A A 2 ) 
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+ 32a 4 A|(Af + 16A A 2 A 2 + B6X 2 X 2 2 ) + 4a 10 A 3 l A 2 (A 4 + 25A A 2 A 2 + 112A§Al) 

+ 8a 8 X 2 Xl(3\j + 40A A 2 A 2 + U0X 2 X 2 2 ) + 4a 6 A 2 ,(A? + 12A A 4 A 2 + U0X 2 XjX 2 2 

+ 448A 3 1 A 3 ,))A 3 s? - (a 2 A + 2A 2 ) 3 (a 16 A 8 1 + 256A 8 , + 16a 12 A^A 2 (A? + 7A A 2 ) 

+ a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A 2 3 (3A 2 + 16A A 2 ) + 32a 4 A^A 4 + 16A A 2 A 2 + 56X 2 Xl) 

+ 4a 10 A 3 ,A 2 (A 4 + 25A A 2 A 2 + 112X 2 X 2 2 ) + 8a 8 A^A^(3A 4 + 40A A 2 A 2 + U0X 2 X 2 2 ) 

+ 4a 6 A 2 ,(Af + 12A A 4 x A 2 + 140A§A?A| + 448A 3 ,A 3 ,))A 3 sis 2 , - 2aAiA 2 (a 4 A 3 1 

+ 4a 2 X 2 X 2 + 2A 2 (A 2 + 2A A 2 ))(a 4 A 2 1 + 4A 2 , + a 2 (A 2 + 4AoA 2 ))(a 10 A^ + 64A^ 

+ a 8 Ag(A 2 + 12A A 2 ) + 8a 2 A^(A 2 + 18A A 2 ) + 2a 6 A 2 1 A 2 (3A 2 + 28A A 2 ) 

+ 2a 4 A 2 (A 4 + 6A A 2 A 2 + 64A 2 1 A 2 ,))s 3 - 2aAi(a 2 A + 2A 2 ) 2 (a 16 A 8 1 + 256A^ 

+ 16a 12 A^A 2 (A 2 + 7A A 2 ) + a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A|(3A 2 + 16A A 2 ) 

+ 32a 4 A^(A 4 + 16A A 2 A 2 + 56AqA 2 ) + 4a 10 A 3 ,A 2 (A 4 + 25A A 2 A 2 + 112X 2 X 2 2 ) 

+ 8a 8 AgA^(3A 4 + 40A A 2 A 2 + 140A^A|) + 4a 6 A|(A? + 12A A 4 A 2 + 140A§A?Al 

+ 448A 3 A 3 ))A 3 s 2 s 3 - a 2 A 2 (a 2 A + 2A 2 )(a 16 A 8 + 1280A 8 + a u X 6 (X 2 + 20A A 2 ) 

+ 8a 12 AfjA 2 (3A 2 + 22A A 2 ) + 64a 2 A^A 2 + 68A A 2 ) + 32a 4 A^(3A 4 + 44A A 2 A 2 

+ 200A 2 1 A 2 ,) + 4a 10 A 3 1 A 2 (2A 4 + 47A A 2 A 2 + 220A 2 1 A 2 ,) + 8a 8 A§A|(7Af + 88A A 2 A 2 

+ 340A 2 1 A 2 ,) + 4a 6 A|(3A? + 32A A 4 A 2 + 348A 2 ,A 2 A 2 . + 1328A^Ai))A 3 sis| 

+ 4aAiA 2 (a 4 Ag + 4A 2 , + a 2 {X\ + 4A A 2 ))(128A^ - 64a 2 A^A 2 - 6A A 2 ) 

+ a 12 A^(2A A 2 - A 2 ) - 8a 4 X 3 2 {Xj + 18A A 2 A 2 - 60AqA 2 ) + a 10 A{j(24A^A| - X\ 

- 12A A 2 A 2 ) + 2a 8 A 2 1 A 2 (60A 2 1 A 2 , - 3A 4 - 28A A 2 A 2 ) - 2a 6 A 2 (A 2 - 2A A 2 ) 
x (A 4 + 8A A 2 A 2 + 80A 2 A 2 ))A 3 s 3 - a 3 Ai(a 2 A + 2A 2 ) 2 (a 4 A 2 + 4A 2 

+ a 2 (A 2 + 4A A 2 ))(a 10 A^ + UX\ + a 8 A|](A 2 + 12A A 2 ) + 8a 2 A|(A 2 + 18A A 2 ) 
+ 2a 6 A 2 ,A 2 (3A 2 + 28A A 2 ) + 2a 4 A 2 (A 4 + 6A A 2 A 2 + 64A 2 A 2 0)A 2 s 2 s 3 

- 2a 4 A 2 (a 2 A + 2A 2 )(a 4 A 2 1 + 4A 2 , + a 2 (A 2 + 4A A 2 ))(a 10 A^ + 96A| 

+ 8a 6 AgA 2 (A 2 + 9A A 2 ) + 16a 2 A^A 2 + 13A A 2 ) + a 8 A|](A 2 + 14A A 2 ) 

+ 2a 4 A 2 (A 4 + 10A A 2 A 2 + 88A^Ai))A|sis 4 + a 3 Ai(a 4 A§ + 4A 2 , + a 2 (A 2 + 4A A 2 )) 

x (128A/J - 96a 2 A|(A? - 4A A 2 ) + a 12 A^(2A A 2 - A 2 ) + a 10 Xl{24X 2 X 2 2 - X\ 

- 14A A 2 A 2 ) + 2a 8 A 2 (60A 2 ,A 2 , - 5A 4 - 36A A 2 A 2 ) + 8a 4 A|(60AgA^ - 3A 4 

- 26A A 2 A 2 ) - 2a 6 A 2 (A? + 14A A 4 A 2 + 88A^A 2 A^ - 160A|]A|))A|s| 

- a 5 Ai(a 2 A + 2A 2 ) 4 a 4 A 2 1 + 4A 2 , + a 2 {X\ + 4A A 2 )) 2 x {\ls\s\ - 2a 6 A 2 (a 2 A + 2A 2 ) 3 
x (a 4 A 2 , + 4A 2 , + a 2 (A 2 + 4A A 2 )) 2 A 3 ; si^ - a 7 A?(a 2 A + 2A 2 ) 2 (a 4 A 2 1 + 4A 2 , 

+ a 2 (A 2 + 4A A 2 )) 2 A 3 4 = 0, 
-2A 2 (a 2 A + 2A 2 ) 3 (a 16 A 8 1 + 256A 8 , + 16a 12 A^A 2 (A 2 + 7A A 2 ) + a 14 A(j(A 2 + 16A A 2 ) 
+ 64a 2 A^(3A 2 + 16A A 2 ) + 32a 4 Af (A 4 + 16A A 2 A 2 + 56A§A|) + 4a 10 A 3 ,A 2 (A 4 
+ 25A A 2 A 2 + 112X 2 X 2 2 ) + 8a 8 A§A|(3Af + 40A A 2 A 2 + MOA^A 2 ,) + 4a 6 Al(A? 
+ 12A A 4 A 2 + 140A§Af A| + 448A 3 ) A 3 >))n - 2AiA 2 (a 2 A + 2A 2 ) 2 (a 16 A 8 l + 256A 8 , 
+ 16a 12 A^A 2 (A 2 + 7A A 2 ) + a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A|(3A 2 + 16A A 2 ) 
+ 32a 4 A|(A 4 + 16A A 2 A 2 + 56A^A^) + 4a 10 A 3 ,A 2 (A 4 + 25A A 2 A 2 + 112X 2 X 2 2 ) 
+ 8a 8 AgA^(3A 4 + 40A A 2 A 2 + 140A^A|) + 4a 6 A|(A? + 12A A 4 A 2 + UOA^A 2 , 
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+ 448AgA|))si + 2a\ 2 (a 2 X + 2A 2 )(512A A^ + a 16 A§(2A A 2 - Xf) - 128a 2 A^(A 4 
+ A A 2 A 2 - WX 2 X 2 2 ) + a u X 6 (32X 2 X 2 2 - A 4 - 16A A 2 A 2 ) + 2a 12 A[jA 2 (112A 2 1 A 2 , 

- 9Xf - 52A A 2 A 2 ) - 32a 4 A|(A? + 13A A 4 A 2 + 16A§Af Al - 112AgA|) 

+ 4a 10 A 3 1 A 2 (224A 3 1 A 3 : - 90AgA 2 A^ - A? - 28A A 4 A 2 ) + 4a 8 A 2 1 A 2 ,(560A 3 ,A 3 ; - 7Af 

- 84A A 4 A 2 - 180A^A 2 Al) - 4a 6 A 2 ,(A? + 14A A?A 2 + 132A 2 , x A 4 A 2 , + 208A 3 1 A 2 A 2 ; 

- 896AgA£))«3 - «(a 2 A + 2A 2 ) 3 (a 16 A 8 + 256A 8 , + 16a 12 A^A 2 (A 2 + 7A A 2 ) 

+ a 14 A(j(A 2 + 16A A 2 ) + 64a 2 A^(3A 2 + 16A A 2 ) + 32a 4 A 2 (A 4 + 16A A 2 A 2 + 56AqA 2 ) 
+ 4a 10 A 3 1 A 2 (A 4 + 25A x A 2 A 2 + 112A^A 2 ) + 8a 8 X 2 X 2 2 (3Xf + 40A A 2 A 2 + 140AgA 2 ) 
+ 4a 6 Al.(A? + 12A A 4 A 2 + U0X 2 XjX 2 2 + 448A 3 1 A 3 ; ))A3S 2 S3 - 2a 4 Ai(a 2 A + 2A 2 ) 2 
x (a 14 A 8 + 128A^(A 2 + 2A A 2 ) + a 12 A(j(A 2 + 16A A 2 ) + 4a 8 A 3 A 2 (A 2 + 5A A 2 ) 
x (A 2 + 20A A 2 ) + 4a 10 A^A 2 (4A 2 + 27A A 2 ) + 32a 2 A|(A 4 + 13A A 2 A 2 + 26A 2 ,A 2 .) 
+ 4a 6 A 2 1 A 2 .(7A 4 + 78A A 2 A 2 + 220A 2 ,A 2 ,) + 4a 4 A^(Af + 14A A 4 A 2 + 128A^A 2 A^ 
+ 288A 3 1 A 2 ! ))A 3 sis 2 ; + a(512a 2 A^(A 2 - 16A A 2 ) - 2048A^1 + a 20 A^(2A A 2 - A 2 ) 
+ a 18 A 7 (2A A 2 - A 2 )(A 2 + 16A A 2 ) + 256a 4 A/^A 4 + 9A A 2 A 2 - 5AX 2 Q X 2 2 ) 
+ 8a m X G X 2 {27X 2 X 2 2 - 2Xf - 9A A 2 A 2 ) + 32a 6 A^(Af + 14A A 4 A 2 + 128A^A 2 A| 

- 384AgA|) + 16a 12 AqA 2 x (84A^A| - 2A? - 15A A 4 A 2 + 16AgA?A|) - 32a 8 A A^(A? 

- 4A A 4 A 2 - 116A^A 2 A^ + 168A 3 A|) + 4a 14 A 4 ;A 2 (192A 3 1 A 2 ! — Af — 23A A 4 A 2 

- 32AqA 2 A 2 ) + 4a 10 A A 2 A 2 ,(432A 3 1 A 2 ! - A? - 18A A 4 A 2 - 60AjjA?Al))A 3 siJ 

+ 2a 3 A 2 (a 2 A + 2A 2 ) 4 (a 4 A 2 1 + 4A 2 . + a 2 (A 2 + 4A A 2 ))(a 6 A 3 1 + 12a 2 A A 2 + 8Xl 

+ a 4 A (6A A 2 - A 2 ))A 2 s 2 s^ + 8a 4 A!A 2 (a 2 A + 2A 2 )(a 4 A 2 1 + 4A| + a 2 (A 2 + 4A A 2 )) 

x (a 10 A^ + 10a 8 A^A 2 + 32A^ + 8a 2 Xl{X\ + 10A A 2 ) + 2a 6 A 2 ,A 2 (A 2 + 20A A 2 ) 

+ a 4 A 2 (A 4 + 8A A 2 A 2 + 80A 2 Ai))A|s lS | - 2a 3 A 2 (a 4 A 2 l + 4A 2 . + a 2 (A 2 + 4A A 2 )) 

x (a 14 A£ + 128X 7 2 - 32a 2 A^A 2 - 14A A 2 ) + 2a 12 A^(7A A 2 - Xj) - 24a 4 A^(A 4 

+ 4A A 2 A 2 - 28AqA 2 ) + a 10 Ag(84A^A| - A 4 - 18A A 2 A 2 ) + 2a 8 A 2 1 A 2 (140A 2 ,A 2 , - 5A 4 

- 32A A 2 A 2 ) - 4a 6 A 2 (A? + 7A A 4 A 2 + 28A 2 ,A 2 A|. - 140A 3 A|))A^1 + a 5 (a 2 A + 2A 2 ) 
x (a 4 A 2 , + 4A 2 . + a 2 (A 2 + 4A A 2 )) 2 (a 6 A 3 1 + 6a 4 X 2 X 2 + 8A 3 , 

+ 2a 2 A 2 (A 2 + 6A A 2 ))A 3 s 2 s| + 2a 6 Ai(a 2 A + 2A 2 )(a 4 A 2 + 4A 2 + a 2 (A 2 + 4A A 2 )) 2 
x (a 6 A 3 , + 6a 4 X 2 X 2 + 8A| + 2a 2 A 2 (A 2 + 6A A 2 ))Aisis^ + a 7 Xj{a 4 X 2 + 4Al 
+ a 2 (A 2 + 4A A 2 )) 2 (a 6 A 3 1 + 6a 4 A 2 A 2 + 8A| + 2a 2 A 2 (A 2 + 6A A 2 ))A 3 4 = 0. 

D Appendix 

The stationary solutions of the system (|16|) : 



'2a 2 A 2 + A 2 - X lv / Aa 2 X 2 + X{ (Ai + V 4 " 2 ^ + A i) \Ao - hs 2 2 

si 



2v / 2a 2 Agv / As 



'2a 2 A 2 + A 2 - Ai V / 4a 2 A 2 + A 2 v / Ao-A 3 s 2 

S3 = 0, r\ = — = — . 

2V2a 2 A 0v ^ 



(Ax - V4a 2 A 2 + X()s 2 1 
r2 = 2^X ' r3 = ° ' 
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si 



2a 2 A 2 + A 2 - Ai v / 4a 2 Ag + A^(Ai + V4a 2 A 2 + A 2 ) V^o - A 3 s 2 



2 V / 2a 2 A 2 v / A 3 ~ 

A /2a 2 A 2 + A 2 -A : 

S3 = 0, ri = 



2a 2 A 2 , + A 2 - AiV4a 2 Ag + A=fVAo-A 3 s^ 



2 v / 2a 2 A ^ 



(Aj - y^A 2 + A 2 )s 2 
r2 = 2^X ' rs 



Here A 2 = -(2a 2 A 2 + A 2 + A lx /4a 2 A 2 + A 2 )/(4A C 



' 2a 2 A 2 + A 2 + Ai V4a 2 A 2 + A 2 (Ai - y / 4a 2 A^ + Af) ^A - A 3 s 2 
Sl ~ 2V2a 2 A 2 VAs" : 



/2a 2 A 2 + A 2 + Ai V / 4a 2 A 2 + A 2 v / Ao-A 3 s 2 

S3 = o, n 



2^2a 2 A 



o 1 



(A 1 + x /4a 2 A 2 + A 2 )s 2 



r2 = 2^ ' r3 = ° 



2a 2 A 2 + A 2 + Ai v / 4a 2 A 2 + A 2 (Ai - v / 4a 2 A 2 + A 2 ) ^An - A 3 s 2 

si 



2^2a 2 A 2 VA 

A /2a 2 A 2 + A 2 + A : 

«3 = o, n = 



2a 2 A 2 + A 2 + Ai y/4a 2 >$ + X{y/\o- A 3 s 2 



2^2a 2 A VA 



3 



(Ai + v / 4a 2 A 2 + A 2 )s 2 
2a 2 An 



r 2 = ^ , ^3 



Here A 2 = -(2a 2 A 2 + A 2 - X 1 v / 4a 2 A 2 + A 2 )/(4A ). 
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